Lubrication approximation for micro-particles moving along parallel walls 
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Lubrication expressions for the friction coefficients of a spherical particle moving in a fluid between 
and along two parallel solid walls are explicitly evaluated in the low-Reynolds-number regime. They 
are used to determine lubrication expression for the particle free motion under an ambient Poiseuille 
flow. The range of validity and the accuracy of the lubrication approximation is determined by 
comparing with the corresponding results of the accurate multipole procedure. The results are 
applicable for thin, wide and long microchannels, or quasi-two-dimensional systems. 
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Microfluidic devices are commonly used for particle ma- 
nipulation and separation, such as biological cell sorting, 
on-chip hydrodynamic chromatography 111], electrophore- 
sis 0, and many other applications [1, 13, H, S 0, B B 
[Tol fllj . Particles often move along microchannels whose 
smallest dimension is comparable with the particles' size, 
and the other ones are much larger [12, BI|. For such a 
quasi-two-dimensional system, it is of interest to theoret- 
ically determine the particle velocity, solving the Stokes 
equations for the fluid motion between two parallel infinite 
walls. 

Particle velocity can be evaluated numerically, by the 
boundary-integral method [3, or the multipole expan- 
sion However, as pointed out in Ref. [Hl |. near-contact 
explicit asymptotic expressions would be useful, in analogy 
to the widely applied lubrication expressions for pair s of 
very close solid spherical particles of different sizes [17| , in- 
cluding the limiting case of a sphere close to a single plane 

wall [SB!. 

Therefore the goal of this paper is to derive lubrication 
expressions for a spherical particle moving along two close 
parallel solid walls, and check their accuracy. 

A spherical particle of radius a moving in a Stokes flow 
between two parallel solid walls is considered, as illustrated 
in Fig. [TJ Distances are normalized by the particle radius 
a. The wall at z = is labeled 1, and the wall at z = h > 
is labeled 2. The sizes of the gaps between the walls and 
the particle surface are small, t\ <C 1 and £2 <^ 1- 

The fluid velocity v and pressure p satisfy the Stokes 
equations [2fj|. [2l|. 



Vp = 0, V 



0. 



(1) 



with the stick boundary conditions at the walls and the 
particle surface. 

In the friction problem, there is no ambient flow and the 
particle translates along the channel with velocity (Jji and 
rotates with angular velocity fl y y. The goal it to find the 
hydrodynamic force T x yi and torque C y y exerted by the 
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FIG. 1: The system and notation. 
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particle on the fluid [201 ] . 
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Here £ denotes the dimensionless 2x2 friction matrix, 



c = 



with the friction factors /* a 
found. 

Free motion in Poiseuille flow means that a spherical 
particle is entrained by the ambient Poiseuille flow, 




Az(z-h)/h 2 x. 
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In the absence of external forces, the particle translates 
along x and rotates along y. The translational and angular 
velocities, normalized by u ma x and v max /a, respectively, 
have the form, 
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with the mobility matrix obtained from Eq. ([3]), 

M = C\ (6) 
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and the dimensionless friction factors f£ x , dP yx determined 
by the force J-Jx and torque C^y exerted by the Poiseuille 
flow on immobile sphere [2J] in the following way, 



known 18, 19, 22, 29], 
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In this work, the functional dependence of the friction 
factors and the particle's free-motion velocities on t\ and 
e 2 has been determined numerically. The accurate theoret- 
ical method [25|, [26[ of solving the Stokes equations, based 
on the multipole expansion [271 ] ■ has been applied. For a 
system of particles between two parallel walls, this method 
involves expanding the fluid velocity field into spherical and 
Cartesian fundamental sets of Stokes flows. The spherical 
set is used to describe the interaction of the fluid with 
the particles and the Cartesian set to describe the interac- 
tion with the walls. At the core of the method are trans- 
formation relations between the spherical and Cartesian 
fundamental sets. The transformation formulas are used 
to derive a system of linear equations for the force multi- 
poles induced on the particle surfaces. The coefficients in 
these equations are given in terms of lateral Fourier inte- 
grals corresponding to the directions parallel to the walls. 
These equations are truncated at an multipole order L [53] 
and solved numerically, using the algorithm and the FOR- 
TRAN code described in Refs. [TJ, [M [26| and available at 

To derive the lubrication approximation, the multipole 
algorithm described above has been first used to compute 
the friction matrix C, given by Eq. (|3|). The friction matrix 
has been evaluated as the sum of two terms, 



fLx(^) = - In £; + 0.95429 -^-e i hxe i + 0.42945 e*, 
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c r lyy (e t ) = --lne* + 0.37089- —e^ne, + 0.34008e 



The task is to derive near-contact approximation of the 
remaining two- wall contribution £ 12 , regular in e^. The 
components of Ci2( £ ii e 2) are linear functions of t\ + e 2 
or ei — £2, depending on symmetry. To determine them 
explicitly, values at the contact, Ci2(0j0), have been first 
evaluated for all the subsequent multi pole orders L < 404, 
and extrapolated to L = oo, see Refs. [24l. W\ for the basic 
concepts. Next, the multipole numerical codes with a very 
high multipole order 402 < L < 404 have been applied 
to evaluate Ci2( e i; e 2) as a function of < e\ < 0.02 for 
e 2 = 0. The following linear relations have been found, 



/i2,xx(ei=e2) « 0.52805- 0.64 (ei + e 2 ), (15) 

(ei,e 2 ) « 0.175 (ei - e 2 ), (16) 
(ei,e 2 ) ~ -0.06643 + 0.08 (ei + e 2 )- (17) 



"12, yx 



with the accuracy 0.52805 ± 0.0001, -0.64 ± 0.01, 0.175 ± 
0.005, -0.06643 ± 0.00002 and 0.08 ± 0.01. 

Finally, the two- wall lubrication asymptotics of the total 
friction matrix ([3]) is the following, 



C(ei,e 2 ) = Cite) -Co. +Ci 2 (ei,e 2 ). (8) f t^ lA3m _ 8 _ 64 + 



The first one, the single-wall superposition, is the sum of 
the one-particle friction matrices in a fluid bounded by the 
wall i only, 
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minus the one-particle friction matrix in unbounded fluid, 

(10) 
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The second one, Ci 2 > is the two- wall contribution, 

Cl 2 = ( I 1 ?'™ 'f 12 ^ ) ■ (11) 

The single- wall near-contact expressions for Ci( e i) are 
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C U~ _ io ln ^ ~ 25l)^ lln 61-62 ln£2 ) + - 276 ( £ i- £ 2), 
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c r w « -0.32465 - - ln( ei e 2 ) - — (e a lnei + e 2 lne 2 ) 
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+ 0.42( £l + e 2 ). 
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For applications, the most important is the translational 
coefficient f* x . For the channel widths h = 2.1 — 2.5, the 
logarithmic term, — ^log(eie 2 ), give only 68.5-50.5% of 
the total lubrication expression (fT5|) . The constant term, 
1.4366, gives the additional 30.8-49.1% (the other terms are 
practically negigible). Note that to evaluate this constant 
precisely, it is essential to compute the two-wall contribu- 
tion (fTS")) at the contact, which is as large as 0.52805, i.e. 
10.0-7.1% of the total lubrication expression lfl8]). 

Derivation of the lubrication formulas for the free mo- 
tion in Poiseuille flow will now be discussed. First, the 
near-contact approximations to the force / x ' x (ei,e 2 ) and 
torque c v yx (t\, e 2 ) exerted by the ambient Poiseuille flow on 
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the motionless particle are constructed, using the accurate 
numerical results. In the absence of relative motion, there 
is no non-analytic terms. For ei = £2 = 0, the force is eval- 
uated for subsequent multipole orders L < 404. The linear 
relation of the force on L~ 3 , given in Ref. (24|, is now used 
to extrapolate the results until L = 00, and to determine 
f£ x (0,Q). The next term is linear in (ej + e 2 ). To de- 
termine the coefficient, the function /f x (ei,0) is evaluated 
with 402 < L < 404, and the linear dependence is found in 
the range < t\ < 0.005. Similar procedure is repeated for 
cP rj ,(ei,0). In this case, however, quadratic terms, propor- 
tional to (e\ — el), are needed in addition to those linear in 
(ei — e 2 ), because, owing to symmetry, €^(0,0) = 0. The 
near-contact asymptotics reads, 

fl x » 2.67817756- 1.1035( ei + e 2 ), (21) 
(? yx « -0.541( ei -e 2 )+0.5(ef-e2). (22) 

Finally, lubrication approximation of the free motion ve- 
locities is obtained from Eqs. (O-©, with the use of the 
near-contact expressions (fT8" |) - (j2"2")l . 

The range of validity of the above lubrication approxi- 
mation has been determined by comparison with the nu- 
merical results. First, motion in a narrow channel of a fixed 
width, h = ei + e 2 + 2, with h = 2.1, 2.2, 2.3, 2.4, 2.5, has 
been examined. The errors of the near-contact expressions 
are listed in Table fl] The absolute errors of f xx , fg x and u 

TABLE I: The maximal absolute and relative errors of the 
lubrication approximation for different channel widths h. 



h 


2.1 


2.2 


2.3 


2.4 


2.5 


^ fxx 


0.004 


0.02 


0.04 


0.06 


0.09 


5c* 


0.0003 


0.002 


0.005 


0.008 


0.01 


8c r 


0.007 


0.02 


0.045 


0.07 


0.10 


XtP 

u J XX 


0.006 


0.02 


0.05 


0.085 


0.13 




0.0003 


0.002 


0.007 


0.015 


0.03 


5u 


0.0008 


0.003 


0.007 


0.013 


0.02 




0.00007 


0.0007 


0.003 


0.007 


0.01 


^ fx x / fx X 


0.09% 


0.4% 


1% 


2% 


3% 


^Cy x j f XX 


0.005% 


0.04% 


0.1% 


0.2% 


0.3% 


Sc r lc r 


0.2% 


0.9% 


2% 


3% 


5% 


gfP try 

w J XX 1 J XX 


0.2% 


1% 


2% 


4% 


5.5% 


8c p IF 


0.01% 


0.1% 


0.3% 


0.7% 


1.3% 


Su/u 


0.15% 


0.5% 


1% 


2% 


3% 


5uj/u 


0.01% 


0.1% 


0.5% 


1% 


2% 



are the largest in the center of the channel, of (f x and c yy - 
at the contact, and of cP xx and w - in between. Note that the 
physical relative errors are Scy X /f xx , 8c^ x / 'fg x and dco/u 
rather than 5c yx /c yx , Sc^ x /c^ x and Su>/lj, respectively. 

In Fig. [2j the accurate friction factors and free motion 
velocities are plotted as functions of ei, together with their 
lubrication approximation, given by Eqs. (fl"8)) - (f2"2"]) and ([5])- 
((6]). Even for a relatively wide channel with h = 2.3, the 
errors of the near-contact physical expressions do not ex- 
ceed 2%. 
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FIG. 2: Friction coefficients (top) and translational and angular 
velocities of a freely moving sphere (bottom) in the channel of 
width h — 2.3. Solid lines: accurate; dashed line: near-contact. 

In Fig. [3l the maximal free motion velocities, u m — 
maxw(ei) and u) m = maxw(ei) as well as the position of 

the largest rotation, e m : u>(e m )=u> m , are plotted versus the 
channel width h. Note that u is the largest at the channel 
center, and u - very close to the wall. For h < 2.5, the lu- 
brication approximation of the maximal particle-velocities 
is still reasonably accurate, with 5u m /u m , Stu m /u m < 3%. 
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FIG. 3: Maximal translational and angular velocities of a freely 
moving sphere, u m and uj m , and the position of the maximal 
rotation, e m , versus the channel width h. Solid lines: accurate; 
dashed line: near-contact. 

Finally, velocity profiles of a freely moving sphere are 
compared for different channel widths. In the computation, 
L = 40, with the 10" 4 absolute accuracy. Both trans- 
lational and rotational velocities of the sphere, u and u), 
are normalized by u m , and the gap e\ between the sphere 
surface and the lower wall - by the sum of both gaps, 
e = £! + e 2 - The results are plotted in Fig.0J The narrower 
channel, the closer to the plug flow translational velocity 



is observed, and the smaller ratio of u>/u. 

Concluding, in this paper the near-contact approxima- 
tion for a solid-sphere motion along parallel hard walls has 
been derived, with at least 2% precision for the distance 
between the walls up to 2.3 radii. 
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FIG. 4: Translational and rotational particle-velocity profiles 
for channels widths h = 2 + e, with e = 0.1,0.2,0.3,0.4 and 0.5 
(top-down for u/u m and down-top for u/u m ). 
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